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Abstract. We propose a new approach to the combinatorial interpretations of linearization 
coefficient problem of orthogonal polynomials. We first establish a difference system and then 
solve it combinatorially and analytically using the method of separation of variables. We 
illustrate our approach by applying it to determine the number of perfect matchings, de- 
rangements, and other weighted permutation problems. The separation of variables technique 
naturally leads to integral representations of combinatorial numbers where the integrand con- 
tains a product of one or more types of orthogonal polynomials. This also establishes the 
positivity of such integrals. 
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1. Introduction 

In the late 1960's Askey formulated several conjectures about the nonnegativity of integrals 
of products of orthogonal polynomials times certain functions. An excellent survey of the 
research in this area, which was spearheaded by Askey, is Askey's CBMS lecture notes [3, see 
also [T]. In the 1970's it was realized that some of the integrals considered by Askey and his 
coauthors have combinatorial interpretations. Even and Gillis [11] showed that the number of 
derangements of sets of sizes ni, n2, . . . , Um is 

/"OO ^ 

Jo 

where L„(3;)'s are the simple Laguerre polynomials, while Azor, Gillis, and Victor [7] and inde- 
pendently Godsil [16] showed that the number of perfect matchings of sets of sizes ni, re2, . . . , rim 
is 

_^,2 m 

2-(ni+-+n™)/2 / ^TJh (x)dx, 

Jr ^ fj^ 

where ff.„(a;)'s are the Hermite polynomials. Askey and Ismail [4] used the MacMahon Master 
theorem to give a systematic combinatorial treatment of the integrals of products of the classical 
polynomials with respect to certain measures. One of them generalized the Even and Gillis 
result to Meixner polynomials. Foata and Zeilberger [12] considered the general Laguerre 
numbers 

Jo r(a + 1) fj- ' 
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where L}f (x)'s are the Laguerre polynomials. Zeng and, Kim and Zeng extended this study 
to all Sheffer polynomials in [Ml [Ml [35] . 

In their combinatorial study of integrals of products of orthogonal polynomials Askey and 
Ismail [1] pointed out another source of positivity results. Recall that a system {Qn{x)} of 
birth and death process polynomials, [20], [T71 §5.2], is generated by 

Qo{x) = l, Qi{x) = [bo + do - x]/do, 

-xQn{x) = bnQn+l{x) + dnQn~l{x) - {bn + dn)Qn{x), 

where {bn} and {dn} are the birth and death rates, respectively, are such that 

(1.3) A„, > 0, n > 0, and (i„ > 0, n > 0, do > 0. 

Karlin and McGregor [2^ showed that the probability to go from state (population) m to state 
(population) n in time t is given by 

(1.4) pm,n{t) = " " " r e-"'*Q™(x)Q„(x) d^(x), t > 0, 

"1"2 • • • "n Jo 

where fi is the orthogonality measure of {Qnix)}- This proves that 

POO 

(1.5) / e-^*Q^(x)Q„(x)(i/i(x) >0. 

Jo 

The Laguerre polynomials correspond to 6„ = n + 1, d„ = n + q. Thus 

/•oo 

(1.6) / e-^*2;"e-^L(^)(2;)L(f)(x)(ix > 0, a > 0, t > 0. 
Jo 

This and the derangement number (II. ip motivated us to consider the combinatorial interpre- 
tation of the numbers 

(1.7) ^W(m,n,.) = ^7 / ^Ll^\^)LtHx)x'^e-^dx. 

T{a + 1) Jo s\ 

One important tool used in the combinatorial study of the integrals of orthogonal polynomi- 
als is MacMahon's Master theorem and its /3-extension due to Foata-Zeilberger [12j . When the 
;0-extension of MacMahon's Master theorem is combined with the exponential formula |30[l33j. 
all the known combinatorial interpretations of the linearization coefficients of the orthogonal 
Sheffer polynomials can be deduced by computing their generating functions. Another way to 
gain insight into the combinatorial interpretation of the linearization coefficients is from their 
corresponding moment sequences, see [M l [32 1 [35]. 

(1.8) M„= / x"(i/i(x). 

Jr 

However the generating function approach fails when one tries to extend the previous results to 
their q-analogues, even though a conjecture for the combinatorial interpretation is formulated. 
For example, an important g-analogue for the linearization coefficients of Hermite polynomials 
was given by Ismail, Stanton and Viennot [18], but their proof remains difficult. We are grateful 
to a referee for pointing out that Effros and Popa rediscovered the Ismail-Stanton-Viennot 
result in [10]. Another proof due to Anshelevich [2] uses stochastic processes, and is also far 
from being elementary. Our paper provides a fresh approach to linearization questions. Indeed, 
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one of the main results of this paper is to give an elementary proof of the Ismail-Stanton- 
Viennot result. 

Separation of variables is a standard technique to solve linear partial differential equations. 
The idea is to seek solutions which are products of single variables then by the principle of 
linear superposition the general solution is a linear combination of these products. The only 
problem left is to use initial and boundary conditions to determine the coefficients. This 
technique can be used to solve difference or differential equations. One important application 
of this method is to solve the Chapman-Kolmogorov equations for birth and death processes, 
see \n\ §5.2]. The latter equations is a system of differential equations in time and partial 
difference equations in two discrete variables whose solution is given by (jl.4p . 

In this paper we show how the separation of variables gives integral representations for 
solutions of certain combinatorial problems. 

Our approach is explained in detail in Section 2. The integrands in our integral represen- 
tations are constant multiples of products of orthogonal polynomials times a measure with 
respect to which the polynomials are orthogonal. The integral representations arise naturally 
through separation of variables of the solution of systems of difference equations satisfied by 
the combinatorial numbers. We may reverse the process by starting with integrals of prod- 
ucts of orthogonal polynomials times their orthogonality measure and reach the combinatorial 
numbers. Some of these integrals arose in problems involving linearizations of products of 
orthogonal polynomials where the focus of attention was their nonnegativity [3 El]. Most of 
the positivity results originated from work by Askey and his coauthors in the late 1960's and 
1970's. For references we refer the interested reader to Askey's monograph [3], and to Ismail's 
book [17]. 

The integral representations studied in this work are of the form 

„ m 

(1.9) / fix)T\pn,{XjX)dtl{x), 

where /i is a discrete or absolutely continuous measure and / is some integrable function. Ismail 
and Simeonov |19j studied the large k behavior of integrals of the form (II. 9p when the rij's 
are all equal. Since the integral in (II. 9p represents the number of ways a certain configuration 
occurs, one can calculate the probability that such configuration occurs. We shall also study 
integrals of the type (|1.9p where the polynomials pnj (x) come from two different families of 
orthogonal polynomials. The positivity results which we establish are not only new but seem 
to be the first of its type. 

The rest of this paper is organized as follows. As we already mentioned in the above para- 
graph our approach is outlined in Section 2, where we characterize the linearization coefficients 
of orthogonal polynomials as the unique solution of some partial differential equations with 
boundary conditions. Then we apply the results of Section 2 to various combinatorial prob- 
lems in Sections 3-8. More precisely, by solving the corresponding partial difference equations 
combinatorially, we deduce the combinatorial interpretations of Hermite and Charlier polyno- 
mials, Laguerre polynomials, Meixner polynomials, Meixner-Pollaczek polynomials, (7-Hermite 
polynomials, g-Charlier polynomials, and g-Laguerre polynomials, respectively. In each case 
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we start with a combinatorial problem involving multisets, deduce a difference equation for the 
combinatorial numbers involved, then identify the orthogonal polynomials which arise through 
the machinery developed in Section 2. Furthermore, in Section 9, we extend the previous results 
to some more general integrals to include the moments, inverse coefficients and linearization 
coefficients. We also compute the corresponding generating functions for the corresponding 
integrals of Lagurre and Meixner polynomials and deduce their combinatorial interpretation 
by applying MacMahon's Master theorem. In Section 10, we give a further extension of the in- 
tegrals of Laguerre and Meixner polynomials. Finally, in Section 11, we prove the crucial step, 
Lemma [8.11 towards the combinatorial solution of the partial difference equations of q-Charlier 
polynomials. 

We follow the standard notation for shifted factorials, hypergeometric functions and their 
g-analogues as in the books [U [T4l [TT] . The work of Koekoek~Lesky~Swarttouw [27] is also a 
standard reference for formulas involving orthogonal polynomials and their basic analogues. 



2. Separation of variables and linearization coefficients 



Let {pn{x)} be a sequence of orthogonal polynomials 

(2.1) / Pm{x)pn{x)d^{x) = Cn,(5m,n, Co = 1- 



The condition (^o = 1 amounts to normalizing total mass of ;U to be 1. Then the polynomials 
{pn{x)} must satisfy a three term recurrence relation of the form 

(2.2) Pn+l{x) = [AnX + Bn] Pn{x) - C„,p„,_i (x) , ?! > 0, 

and we will always assume po{x) := l,pi{x) = Aqx + Bq. Therefore 

^0 

(2.3) Cn = -7-C1C2 ■ ■ ■ Cn- 

We consider the linearization coefficients in the expansion of Y\Y=i Pnji'^jx) in {pnix)}. 
Equivalently we consider the numbers 

(2.4) /(n) := ^(^ J] Pn^J^^jPuM di^{x), 

where n = (ni, . . . , Um), nj is a nonnegative integer for 1 < j < m. We shall use the following 
notation: 

Ip{n) = /(ni, . . . ,nj_i,nj ± l,nj+i, . . .,nm)- 
Moreover we assume that = 1- It is clear that 

lf{0,...,0,n) = XjCi — dn,i + Bo{l-Xj)6n,o, if n = 0,l, and 

A.I 

(2.5) m-l 

1(0, . . . , 0) = 1, /(n) = if < n. 
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Theorem 2.1. The numbers I{n) satisfy the system of difference equations 

(2.6) ^ti'^) - Ujk{'n)I^{n) = [B^^ - n^fc^nj I{n) - Cn^Ijin) + Ujk{n)CnJ^ (n), 

where Ujk{n) = Vj{n)/vk{n) and Vj{n) = An^Xj. 



Proof. For 1 < t < m, we have by (I2.2p 

It{'^)= I [{An^XtX + Bnt)Pnt{^tx) - CntPnt-l{Xtx)\Wpnr{Kx)dll{x) 
„ m 

= vt{n) I xW pnX\-x)d^JL{x) + BnJ{n) - CnJti'^)- 

7=1 

Specializing the above equation at t = j and t = k immediately leads to (j2.6p . □ 



Observe that in the system (|2.6p we assmne j,k > 1. It is more convenient to write (j2.6p in 
the more symmetric form 



1 



Brii B„ 



(2.7) ^— /+(n) — I+(n)= r^Hn) I^j- (n) + —^J- {n 



Vk{n)' 



_Vj{n) Vk{n) 

We will show that the system ()2.7p describes many combinatorial problems. Prom now on we 
will consider different combinatorial problems and derive a system of equations of the type 
(|2.7p for the combinatorial numbers under consideration. Theorems 12.21 and 12.31 identify the 
combinatorial numbers as integrals of products of orthogonal polynomials. 

Theorem 2.2. One solution to 

(2.8) y+(n) - Ujk{n)y+{n) = [B^^ - Ujk{n)Bn^]y{n) - Cn^yj{n) + Ujk{n)Cn^yi^ (n), 
is given by 

„ m 

(2.9) y{n) = / Hpn.iXjx) dv{x), 
for any measure v having finite moments of all orders. 

Proof. We try the separation of variables y(n) = YYJLi^ji^j)- When we substitute in ()2.8p 
we get 

Fj(n, + 1) Cn,Fj{nj-l) B^^ _ Fk{nk + 1) , C„,Ffc(nfc - 1) 



(2.10) J ' ^ + -'^j-jy-j _ ""J ^ ' _^ 

Vj{n)Fj{nj) Vj{n)Fj{nj) Vj{n) Vk{n)Fk{nk) Vk{n)Fk{nk) Vk{n) 

Thus each side of the above equation is a constant independent of j or k, so we denote the 
constant by x. This leads to the difference equation 

Fj{nj + 1) = {XjAn^x + Bn^)Fj{nj) - Cn^Fj{nj - 1), 

and the F's now depend on x. Comparing with (|2.2p and noting that Fj{—1) = and Fj{0) = 1, 
we see that the above recurrence relation has a solution given by Fj{nj) = PnjiXjx) and by 
the principle of linear superposition the function in (12. 9p is a solution. □ 
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Theorem 2.3. The system of equations (12 .Sp and the boundary conditions 

y^{0,...,0,n) = XjCi^ dn,i + Bo{l- \j) 5nfl, if n = 0, 1, and 

(2.11) m-l 

y(0,...,0) = 1, y(n) = if ^ < n, 

i=i 

have a unique solution which is given by (j2.4|) . 

Proof. We know that the multisequence (j2.4p satisfies the system of equations ()2.8p and bound- 
ary condition ()2.1ip . hence a solution exists. The second boundary condition defines y for n > 
and when the rest are zero. The first boundary condition in (12. lip defines y for n and when 
one other entry = 1 and the rest are zero in a unique way. Letting n = (0, . . . , 1, 0, . . . , n), 
n > in (j2.8p with k = m we evaluate y(0, . . . , 2, 0, . . . , n) and by induction we evaluate 
y{0, . . . , rig, 0, . . . , n), 1 < s < m. Next we use (|2.8p to evaluate y for general rir, n and another 
nonzero n^: if n > + 1, then y{n) with nonzero entries in the positions r,m oi n is zero 
when we have 1 in the position s; if n < + 1, we use (j2.8p with j = s,k = m to evaluate 
y. Thus we can evaluate y inductively when n has three nonzero entries. We continue this 
argument until we reach any desired general n. □ 

Remark 2.4. It is important to note that (j2.8p is satisfied by solutions of the form (j2.9p 
where v is any probability measure with finite moments. It is the boundary conditions (j2.1ip 
that force v to be an orthogonality measure of{pn{x)}. 



An important class of orthogonal polynomials is the class of birth and death process poly- 
nomials. They are generated by (jl.2p . These polynomials have only positive zeros so they are 
orthogonal with respect to a probability measure supported on a subset of [0, oo). The idea of 
separation of variables is also used to solve the differential-difference equations describing this 
model, see §5.2 and Theorem 7.2.1 in [17j . Birth and death processes have many applications 
in applied probability and queueing theory. 



An immediate consequence of Theorem 
{Qn{x)} generated by (|1.2p . 



is the following result for the polynomials 



Theorem 2.5. The system of difference equations 



(2.12) 

and the boundary conditions 



brij dn 



d, 



(2.13) 



yUO, . . . ,0,n) = Aj -i 5n,i + 



Afc 



1 + ? 

bo 



y{n) 



-J-Vj (^) + -^Vk (^)' 



if 



n 



2/(0,...,0) = l, y(n) = if 



0, 1, and 

m—l 

< n, 

i=i 
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have a unique solution which is given by 

/■OO ™ 

(2.14) y{n)= / TT Q„^. (A^x) d^(x), 

where fi is an orthogonality measure for the polynomials {Qnix)} in (ll.2p . 

From combinatorial point of view, sometimes it is easier to establish a different kind of 
difference equations from (j2.6p . Since pi{x) = Aqx + Bq and pn+i{x) = [AnX + Bn]pn{x) — 
CnPn-iix), we have 

^0 

(2.15) XjPi{x)pn^{XjX) = —pn^+l{Xjx) 



+ [XjBo - ^j— -Bn, )Pn,{Xjx) + -^Cn^Pn^-l{XjX). 



Substituting in (12. 4j) yields 

XjAn 



(2.16) /(i,n) = -^/;(n) + 



/H + ^%/r(n). 



Subtracting (|2.16p from itself with j replaced by k, we obtain (|2.6p . For the Laguerre polynomi- 
als, g'-Charlier polynomials and g-Laguerre polynomials, we shall first establish combinatorially 
()2.16p before passing to (12. 6p . Finally we have the following result. 

Theorem 2.6. Let n = {ni, . . . ,nm) with m > 1 and ni, . . . ,nm. > 0. Any sequence I{n) 
satisfying the system ()2.16p is uniquely determined by its special values at n = (1, . . . ,1) and 
the symmetry with respect to the indices ui, . . . ,nm- 

Remark 2.7. Evaluating the special values of the linearization coefficients ain = (l,...,l) 

amounts to computing the moments of the corresponding orthogonal polynomials, while the 
boundary condition (I2.13P is much easier to check and does not need the knowledge of the 
moments, though the latter would be a source of inspiration for the linearization coefficients. 



3. Linearization coefficients of Hermite and Charlier polynomials 



In this section we consider the linearization coefficients of Hermite and Charlier polynomials. 
We start with some combinatorial setup, which will also be used in the later sections. 

3.1. Combinatorial definitions. In the sequel, we denote by TWn, n„ and (3„ the set of 
perfect matchings, of partitions and of permutations, respectively, of [n] := {1, 2, . . . , n}. Recall 
that a perfect matching of [n] is just a set partition of [n] the blocks of which have exactly two 
elements. It is often convenient to represent pictorially set partitions and permutations of [n]. 
We first draw n elements on a line labeled 1, 2, . . . , n in increasing order. Then, the diagram 
of a partition of [n] is obtained by joining successive elements of each block by arcs drawn in 
the upper half-plane. Here, we say that two elements i < j m. the block B are successive, or 
more precisely that j follows i, if there is no element p £ B such that i < p < j. We denote by 
{i,j) the arc whose extremities are i and j. The diagram of a permutation a G &n is obtained 



SEPARATION OF VARIABLES IN COMBINATORICS 



9 



12345678 12345678 



1 




Figure 1. Diagrams of, from left to right, the matching M = 14/2 6/3 7/5 8, 
the partition vr = 1 4/2 3 7/5 8/6 and the permutation a = 83746251 

by drawing an arc i — a{i) above (resp. under) the hne if i < cr{i) (resp. i > cr{i))- Arcs are 
always drawn in a way such that any two arcs cross at most once. 

In what follows, we fix an m-tuple of nonnegative integers n = (rii, . . . ,nm) such that 
n = ni + • • • + Um and partition the n balls {1, . . . , n} into m boxes Si, ... , Sm where Sj = 
{rii + • • • + Uj^i + 1, . . . , ni + • • • + rij}, no = 0, for j = 1, . . . , m. We denote by [n] the set 
{1, . . . , n} with underlying boxes /Si, ... , Sm, and the corresponding sets of matching, partitions 
and permutations by 



A partition vr of [n] is said to be inhomogeneous if each block of vr contains at least two elements 
and no two elements in the same block belong to the same box Si {1 < i < m). Similarly, 
a permutation a of [n] is an inhomogeneous derangement if cr(5j) H = for all i G [m]. 
We let /C(n) (resp., V{n) and D{n)) denote the set of inhomogeneous perfect matchings 
(resp., partitions and derangements) of [n]. Note that a set partition (resp., permutation) is 
inhomogeneous if and only in its diagram, there is no isolated vertex and no arc connecting 
two elements in the same box Sj (1 < j < m). For instance, if n = (2,3,3), then in Figure [T] 
the matching drawn is in /C(n) while the partition and the permutation are not in V{n) and 
P(n) (they have isolated points). Inhomogeneous objects are drawn in Figured 



Figure 2. Diagrams of, from left to right, a partition in V{2, 3, 3) and permu- 
tation in P(2, 3, 3) 

3.2. Hermite polynomials and inhomogeneous matchings. The Hermite polynomials 
{Hn{x)}n>o can be defined by one of the following five equivalent conditions: 

(1) (Coefficients) H^{x) = Eo<2fc<n(-l)' 2.,!(n-2fe)! (2^)""''- 

(2) (Generating function) Yl'^=o ^nix)^ = exp(2xi — t^). 

(3) (Orthogonality relation) J^Hm{x)Hn{x)e~^^ dx = 2"'n\ ^ 6mn- 

(4) (Recurrence relation) 2xHn{x) = Hn+i{x) + 2nHn-i{x), with H^i{x) = 0, Hq{x) = 1. 

(5) (Moments) //2n+i = 0, /i2„, = 1 • 3 • • • (2n - l)/2". 

Let K{n) be the number of inhomogeneous perfect matchings of [n]. 



(3.1) 



M{n) ■.= M 



n(n) := n„ and 6(n) := 6, 





1 2 3 4 5 6 7 8 
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Lemma 3.1. For k,j E [m] and k ^ j the numbers K{n) satisfy 
(3.2) K+{n) - K+{n) = nkK-{n) - njK-{n), 

and the boundary condition (I2.1ip with Xj = 1 for all j and AqCi = Ai. 

Proof. Let r G [m] and set Nr = ni + • • • + n^. For any i ^ r, the number of matchings 
in /C,^(n) in which A'r + 1 is matched with an element in Si is clearly niK~{n). This implies 
that for any r G [m], we have 



i=l 



from which we immediately deduce (j3.2p . The boundary conditions in (|2.1ip are obviously 
satisfied. □ 

Theorem 3.2. The numbers K{n) have the following integral representation 



(3.3) K{n) = 2-("i+-+"-)/2 



Proof. When Xj = 1 for all j, and 



m 

e 



'—j=^ TT HnAx)dx. 



An = 1, Bn = 0, Cn = n for all n, 



by Lemma [3. H the numbers K{n) satisfy (12. Sp and (12. lip . On the other hand, the correspond- 
ing orthogonal polynomials Hn{x) satisfy the recurrence relation 

(3.4) xHn{x) = Hn+i{x) + nHn^i{x), with H^i(x) = 0, Ho{x) = 1. 

Hence, these are the normalized Hermite polynomials Hn{x) = 2^"/2//„(x/\/2) and their 
orthogonality relation is 

[ ~ ~ e-^'/2 

(3.5) / Hm{x)Hn{x) dx = n\ 5mn- 

JR V 27r 

Therefore 

K{n)= / ^^r[2-"/2i/ (x/V2)dx, 

which is equal to (|3.3p . □ 

Note that the exponential formula (see |301 Corollary 5.1.6]) implies that 
,3.) ;.(„)^...i^.exp 

ni,...,nm>0 ^ i<j 

Hence (j3.3p can also be proved from the generating function of Hermite polynomials. 
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3.3. Charlier polynomials and inhomogeneous partitions. The Charlier polynomials 
Cn^ (x) can be defined by one of the following five equivalent conditions: 

(1) (Explicit formula) ct\x) = ELo © (fc)^K-«)""'- 

(2) (Generating function) En=oCt\x)^ = e-'^"'(l + t/;)^. 

(3) (Orthogonality) ci^\x)dr^\x)dip^''\x) = a"n!(5 

mni where tp^"'^ is the step function 
of which the jumps at the points x = 0,1, . . . are ^("^(x) = ^-^r~. 

(4) (Recursion relation) cl^J^^{x) = {x — n — a)Cn\x) — anC^'^\(x). 

(5) (Moments) = X]fc=i S{n, k)a^, where S{n, k) are the Stirling numbers of the second 
kind. 

The number of blocks of a set partition vr is denoted by bl(7r). Consider the enumerative 
polynomial of inhomogeneous partitions 

(3.7) C{n;a)= a"W. 

Lemma 3.3. For k,j G [m] and k ^ j the polynomials C{n;a) satisfy 

(3.8) ^/(^; ~ ^fc"!"-' ^) — ('^fc ~ nj)C{n; a) + ankCj^{n; a) — anjCj{'n; a), 
and the boundary condition ()2.1ip with Xj = 1 for all j and AqCi = Ai. 

Proof. Let Nj = ni + ■ ■ ■ -\- nj. The partitions of Vj~{n) can be divided into three categories: 

• Nj + 1 and one element of Sk form a block of two elements, the corresponding generating 
function is an^C^ {n; a); 

• + \ and one element of belong to a block containing at least one another element, 
the corresponding generating function is Yl-K(^v{n)^''^k — "-fcj(^))a^'^'^\ where nkj{TT) is 
the number of blocks in vr containing both elements of Sj and (clearly ?ifcj (7r) = 

• Nj + I is in a block without any element of S^, let Rkj{n;a) be the corresponding 
generating function. 

Thus we have 

(3.9) ^^/C^; = X] ~ nkj{'rT))a^^'^''^ + ankC^{n; a) + Rk,j{ri; a). 

Exchanging k and j in the latter identity and subtracting the resulting identity from the latter 
identity, we obtain (j3.8p in view of the symmetry relation Rkj{n; a) = Rj k{n; a). This relation 
can be easily proved, for instance by observing that a partition in 'P^(n) can be seen as an 
inhomogeneous partition of the union U ^2 U • • • U Sm with Si = [A^t_i + 1, Ni] for i / j and 
Sj = [Nj^i + 1, Nj] U {x} where x is any object which is not in [n]. □ 

Remark 3.4. We can also argue as follows. Let n* = (ni, . . . , nj, 1, n^+i, . . . , n^) with j G 
[m] . R is fairly easy to show that 

(3.10) C{n*;a) = Cj~{n; a) + anjCj{n; a) + njC{n; a). 
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Subtracting the above identity from (I3.10p with j = k yields immediately (I3.8p . We will use 
this argument for the Laguerre and Meixner polynomials in the next sections. 

We can solve the system (j3.8p by applying the method of separation of variables which 
naturally leads to the Charlier polynomials. 

Theorem 3.5. The polynomials C{n;a) have the following integral representation 

POO 

(3.11) C{n;a)= C^:;> {x) ■ ■ ■ C^-l{x)di^^'^\x). 

Jo 

Proof. Clearly (j2.8p reduces to (j3.8p when \j = 1 for all j, and 

An = 1, Bn = —n — a, Cn = an for all n > 0. 
From Lemma 13.31 and Theorem 12.31 we deduce (|3.1ip . □ 

The above formula was first established by Zeng |34j using the generating function and the 
exponential formula. A different proof was given by Gessel [15j using rook polynomials. 

4. Linearization coefficients of Laguerre polynomials 

The shifted factorials are 
(4.1) (a)o = 1, (a)„ = a(a + 1) • • • (a + n- 1), n > 0. 

The Laguerre polynomials are defined by 

and have the generating function 

-xt 



(4.3) 5^4-)(a;)r = (l-r""'exp(- 



^ t 

n=0 



They satisfy the recurrence relation 

(4.4) (n + l)Li%ix) - (2n + a + 1 - x)L(")(x) + (n + a)4"A(x) = 0, 
and the orthogonality 

(4.5) r ^f^M^Hx)L(^-\x)dx =^-^±^5m,n. 
Jo r(a + l) n! 

The moments are 

(4.6) = ^, , / x"+"e-^dx = (a + 1)„. 

r(a + 1) 7o 
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In this section we shall prove the results of Foata and Zeilberger [1^ about the Laguerre 
polynomials through our method of separation of variables. For vr E we let Fixj vr = 

7r(5'.j) n Si for i £ [m]. For an ?n-tuple A = (Ai, . . . , A^), define 

m 

(4.7) L{n;a,A)= ^ (q + 1)^^^^ J](Ai - 1)1^'"^"^^^^'"'"', 

7re6(n) i=l 

where cyc(7r) is the number of cycles of tt. By definition, for an inhomogeneous permutation 
vr S we have |Fixj7r| = 0. Hence, when A = 1 := (!,...,!) the summands in (|4.7|) 

reduce to (a + iyy^(^^ if vr € D{n) and otherwise. Thus, we have 

(4.8) L{n;a,l)= ^(a + l)^^^^. 

TT£V{n) 

Lemma 4.1. For j,k G [m] such that j ^ k the polynomials L{n;a,A) satisfy 

(4.9) AfcL+(n; a, A) - XjL+{n; a, A) 

= [(2nfc + a + l)Xj - {2nj + a + l)\k] L{n; a, A) 

+ nk{nk + a)XjLl{n;a,A) - nj{nj + a)XkLj {n; a, A). 

Proof. Let ng = Aq = 1 and n* = (uq, ni, . . . , n^yi) with Sq = {1*}. Then XjL{n*; a, A) is the 
generating function of cj € (3(n*) such that ct(1*) ^ 1* and the edge 1* — t- cr(l*) is weighted 
by Xj. We show that 

(4.10) XjL{n*;a, A) = L+(n; a. A) - (a + l){Xj - l)L{n; a, A) 

+ 2njL{n; a, A) + nj{nj + a)Lj{n; a, A). 

To do so, we adjoin the element 1* to Sj. Thus (a + l)(Aj — l)L(n;a, A) is the generating 
function of a G S^(n) such that cr(l*) = 1*. Hence, the difference 

L^(n; a. A) - (a + l)(Aj - l)L{n; a, A) 

is the generating function of cr G (3^(n) such that cr(l*) 7^ 1*, moreover the edge 1* — )• 0"(1*) 
is weighted by Xj — 1 if ct(1*) G 5j and Aj otherwise. To compensate the over counting, we 
should add 

• the generating function of a G (3^(n) such that o"(l*) G Sj and the edge 1* — cr(l*) is 
weighted by 1; 

• the generating function of a G (3^(n) such that a~^{l*) G Sj and the edge a~^{l*) — t- 
1* is weighted by 1. 

For any a G (3^(n), we let a = cr(l*) and b = a~^{r). There are four cases to consider. 

(1) a G Sj and b ^ Sj. We can construct such a permutation a as follows: starting 
from a permutation r G 6(n) and choosing a point ^ G Sj, we define cT(a;) = t{x) if 
X 7^ l*,T-i(0, and cj(l*) = ^, o-(r-i(^)) = 1* As the weight of the edge 1* ^ ^ is 1 
and that of t~^{(^) —t- 1* in cr is equal to that of t~^(^) — )■ ^ in t, the weight of a is 
equal to that of r, hence the generating function is njL{n\ a. A). 
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(2) a ^ Sj and h G Sj. Similar to the above case, the generating function is njL{n; a, A). 

(3) a £ Sj and b £ Sj, but a ^ b. Starting from a £ &{n*) we can define the permutation 
r on [n] \ {a} by r(j) = a{j) for j ^ a,b and r(&) = a{a). Clearly cyc(cr) = cyc(r). 
Inversely, starting from a permutation r £ &J (n) there are rijijij — 1) choices for a 
and b. Thus, the corresponding generating function is nj{nj — l)L~(n;Q, A). 

(4) a = b £ Sj. The generating function is (a + l)njLj (n; a. A). 

Summing up the above four cases we obtain (j4.10p . Now, substituting j by k in (|4.10p yields 

(4.11) XkL{n*;a, A) = L+(n; a. A) - (a + l)(Afc - l)L(n; a. A) 

+ 2nfcL(n; a, A) + nk{nk + a)L'^{n; a, A). 

Multiplying (j4.10p and (j4.1ip by Xk and Aj, respectively, and then subtracting, we obtain the 
identity (gj]). □ 

We need to state some preliminary results before proving the main result of this section. Let 
A and B be two disjoint sets of cardinality m and n, respectively. An injection / from A to AUB 
can be depicted by a graph on AL) B such that there is an edge a; — )• y if and only if f{x) = y. 
Hence the connected components of the graph consists of cycles, i.e., x — )• f{x) —)■•••—)• /'(x) 
with /*(x) £ A and /'(x) = x and paths, i.e., x — )• f{x) — t- • • • — ?• f^{x) with /'(x) G B. Let 
cyc(/) be the number of cycles of /. Then, Foata and Strehl [13] proved 

(4.12) Yl ^'^^^ = {P + n)^. 

f:A-^AVjB injection 

Theorem 4.2. The polynomials L{n;a,A) have the following integral representation 

(4.13) L(n;a,A) = (-l)"^+-+"-^^^-^ x^e"- H 4'^)(A,x)dx. 

Moreover, this formula is equivalent to the special A = 1 case. 

Remark 4.3. Let us first explain what we mean by the equivalence of (|4.13p and its special 
A = 1 case. We first prove that the definition (14. 7p implies that L{n; a, A) is given by the 
integral (|4.13p . By taking A = 1 in (|4.13p . the formula reduces to the special A = 1 case. The 
point is that we shall prove that knowing the equality in (]4.13p for A = 1 proves that the two 
sides of (I4.13P are equal via the use of the well-known formula [171 Theorem 4.6.5].' 

(4.14) 4°'(-) = (°+i)..E (;_',-:\.)/ i°'w. 

The formula (j4.13p was first proved by Even and Gillis [11] for a = and A = 1. Foata and 
Zeilberger [12] proved the general case of ()4.13p by introducing the cycles. 

Proof. Clearly (|2.8p reduces to (|4.9p when A^ = 1, and 

An = l, Bn = -{2n + a + l), C„ = n(n + a) 
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for all n. That is, the orthogonal polynomials are the normalized Laguerre polynomials Pn{x) = 
(— l)"'n!Li"^(x), which satisfy the three-term recurrence relation 

(4.15) xpn{x) = Pn+i{x) + (2n + Q + l)pn{x) + n(n + a)pn-i{x), 

and the orthogonal relation 

„Q „ — X 

Pn{x)pm{x)dx = n\{a + 1) 



r(a + i) 

From Lemma 14.11 and Theorem 12.31 we deduce (j4.13p when Am = 1- To recover the general 
Am 7^ 1 case, we can proceed as follows: let E he a subset of Sm with cardinality rim — k, 
we consider the permutations vr of 6(n) such that Fixm('7r) = E. Any such a permutation 
corresponds to a pair (a, r) such that a is the restriction of vr on E, which is an injection from 
E to Sm, and r is a permutation on S*! U • • • U Sm-i U {Sm \ E) defined by t{x) = it{x) if 
7r(x) ^ E and t{x) = tt''{x) where / is the minimum integer such that tt''{x) ^ E. Clearly, the 
correspondence vr i— >• {a, r) is a bijection and the generating function of such permutations is 

(a + 1 + A:)„^_,,(Am, - 1)"'"-'=A^ L{n*^; a. A*), 

where n'!^ = (ni, . . . , Um-i, k) and A* = (Ai, . . . , Am-i, 1)- Applying the result for A^ = 1 
case we obtain 

L{n; a. A) = ^ i^'j^j (a + 1 + A:)„,„_fc(Am - iT^-'xt, a. A*) 



fc=0 



m. ^ Q, 2; m-1 



fc=0 ^' 

Now, invoking the known formula (j4.14p . we deduce (|4.13p . 

It remains to show the special A = 1 case of (j4.13p implies ()4.13p for general A. As in the 
above argument, instead of operating within the last box, applying the same operation to all 
the boxes and using ()4.13p for A = 1 we obtain 

L(n; a. A) = Yl 11 ( ) + ^ + ^i)".-'^'. " 1)"^"'^ A^'^L(k; a, 1) 

(a + l + fc)„^._fc(l-A,)"^-^A^:^(,)^ ^ 



\j — 1 k — / 

Thus, the general formula (j4.13p follows by applying the multiplication formula ()4.14p . □ 
Remark 4.4. The analogue of ()4.13p for Hermite polynomials \25\ Proposition 5.1] reads 

__^2 m m 
(4.16) 2-("^+-+'^'»)/2 / ^Y\Hn{Xix)dx= y TT(a2 _ l)hom,W J5.h2hom,W^ 

•^'^ ^ i=l ^eM{n)i=l 
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where homj(7r) denotes the number of homogeneous edges in Si for 1 < i < m. When = 1, 
the right-hand side of (I4.16P reduces obviously to the number of inhomogeneous matchings of 
[n], so the formula (I4.16P becomes (I3.3p . As the analogue of (I4.14p for Hermite polynomials 
[ni (4.6.33)] is 

(4.17) H^icx) = Y: klin-LyS ^ - ^')'^""''^n-2.(x), 

a similar proof of (|4.16p from l\3.3\] using (|4.17p can be given. We leave this to the interested 
reader. 



5. Linearization coefficients of Meixner polynomials 

The Meixner polynomials are \n\ [27| 

(5.1) Mnix;/3,c) = if3)n2Fi{-n,-x;/3;l-l/c), 
and satisfy the orthogonality relation 

(5.2) Mmix; /3, c)M„(x; /3, c)^ = jf^^"_''' '^m.n, /3 > 0, < c < 1. 

The Meixner polynomials generalize the Laguerre polynomials in the sense 

lim Mn> 

They have the generating function 



lim Mn{x/{1 -c);a + 1, c) = n\L^{x). 



(5.3) V M^[x- /3, c)- = (1 - i/c)"(l - t)-"- 

n=0 

The notation here is slightly different from [TTJ Chapter 6] . The three-term recurrence relation 
is 

(5.4) -xMn{x;l3,c) = cil-cy^Mn+i{x;l3,c) 

- [c(/3 + n) + n](l - c)-^Mn{x; (3, c) + (1 - c)-\p + n- l)nM„_i(x; (3, c). 
The moments are, see (28 | [32 t [35], 

(5.5) Mn(/3, c) = (1 - c)/^ 5: rc^i§^ = ^-^""^rrf 

fc>o ■ ^ ' 

where wex(7r) is the number of weak excedances of vr, i.e., 

(5.6) wex(7r) = 1 < i < n and i < 7r(i)}|. 

Let vr be a permutation of [n]. We say that vr has an excedance (resp. box-excedance) at 
i £ [n] if i < 7r(z) (resp. z E 5^, ir{i) G S'j and j > A:). Denote by exc(7r) (resp. excb(7r)) 
the number of excedances (resp. box-excedances) of vr. Clearly, if vr is an inhomogeneous 
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derangement, then exc(7r) = exc;,(7r). Consider the generating function of the derangements 
with respect to the numbers of cycles and (box-)excedances: 

(5.7) M(n;/3,c)= ^ ^cycW^cxcW^ 

neV{n) 

Lemma 5.1. For any k,j € [m] such that k ^ j we have 

(5.8) M+(n; /3, c) - M+(n; /3, c) = (c + l)(nfc - n,-)M(n; /3, c) 

+ cnfc(nfc + /3 - l)M^-(n; /3, c) - cnj(nj + /3 - l)Mr(n; ^, c), 
and the boundary condition (I2.1ip wzi/i Aj = 1 /or a// j and AqCi = Ai. 

Proof. Let j E [m] and set n* = (ni, . . . , rij, 1, rij+i, . . . , n^)- We first show that 

(5.9) M{n*;^, c) = M+{n; (3, c) + (c + l)M(n; /?, c) + nj(nj + /? - l)Af- (n; /3, c). 

Let ti = ni + • • • + rij + 1 and for each vr G D{n*), let a := 7r(ti) and b := 7r~^(ti). We partition 
the derangements in D[n*) into five categories: 

(1) a ^ Sj and 6 ^ Sj. These derangements can be easily identified with the derangements 
in Dj'{n), so the corresponding enumerative polynomial is Mj~{n;f3,c). 

(2) a ^ Sj and b G Sj. Define the derangement vr' on [n*] \ {u} by TT'{j) = Tr{j) for j ^ b 
and '/r'(6) = a. Clearly cyc(7r') = cyc(7r) and exc(7r') = exc(7r) — 1. Conversely, starting 
with any derangement vr' of [n*] \ {u}, we can recover a derangement vr G D[n*) by 
choosing any element in Sj as b and breaking the arrow b — t- '/r'(6) into b ^ u and 
u —7- TT'{b), so the corresponding enumerative polynomial is cnjM{n; (3, c). 

(3) a G 5j and 6 ^ Sj. Define the derangement vr' on [n*] \ {u} by '/r'(j) = 7r(j) for j / 6 
and vr'(6) = a. Clearly cyc(7r') = cyc(7r) and exc(7r') = exc(7r). As in the case (2), the 
corresponding enumerative polynomial is njM{n;/3,c). 

(4) a = b and a £ Sj. The corresponding enumerative polynomial is cjiujMj {n; f3,c). 

(5) a £ Sj, b £ Sj, and a ^ b. Define the derangement vr' on [n*]\{a, u} by 7r'(j) = 7r(j) for 
j b and vr'(6) = 7r(a). Clearly cyc(7r') = cyc(7r) and exc(7r') = exc(7r) — 1. Conversely, 
starting with a derangement on [n*]\{u, a}, we can reverse this process by choosing any 
element in Sj \ {a} as b. As there are nj{nj — 1) ways to choose two different elements 
a and b in Sj, the corresponding enumerative polynomial is cnj{nj — 1)M~ (n; /3,c). 

Summarizing the above five cases leads to ()5.9p . Specializing (15 .Op aX j = k and then subtract- 
ing the resulted equation from ()5.9p ends the proof. □ 

Theorem 5.2. We have 

(5.10) M(n; /3, c^') = - c)^ ^JJ M„^. (x; /3, c)^|p. 

a;=Oji=l 

Proof. When Aj = 1 for all j, and 

^„ = l-l/c, Bn = l3 + n + n/c, C„ = n(/3 + n - l)/c for all n > 0, 
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by Lemma EH the polynomials (-l)"i+-+"'»M(n; /?, c"^) satisfy (USD and (fZTT]) . Theo- 
rem [231 implies then (jS.lOp . □ 

This formula was first given by Askey and Ismail |4] when /? = 1, and by Zeng [M] for 
general (3. 



6. Linearization coefficients of Meixner-Pollaczek polynomials 
The Meixner-Pollaczek polynomials Pn{x) := Pn{x;5,r]) can be defined bv[^ [27], 



X arctan ( 



t 



\l + 6t 



(6.1) V Pnix; 5, 7?)- = [(1 + 5tf + exp 
They satisfy the recurrence relation: 

(6.2) Pn+i{x; 5, 77) = (x - (<5 + 2n)r/)P„(x; 6, ry) - n(7y + n - 1)(1 + S^)Pn-iix; 6, ry). 
The orthogonality relation is 

(6.3) ^— - / Pnix)Pmix)w{x) dx = {5^ + l)"n!(r?)„5„,, 



j)cxc{o-)^cyc{cr)^ 



J^w{x)dx 
where w{x) = x{x; 5, rf) is given by 

— ^ — j exp(—x arctan (5). 

Recall that a permutation vr of [n] has a c/rop (resp. box-drop) at i £ [n] if f > 7r(f) (resp. 
i S 5fc, 7r(z) G S'j and j < k). Denote by drop(7r) (resp. drop{,(7r)) the number of drops (resp. 
box-drops) of vr. 

The moments of Meixner-Pollaczek polynomials [35] are 

(6.4) finis, ii) = , ] I x^'wix) dx=y^{5 + z)'^'-°pM {5 

where = —1. 

Consider the enumerative polynomial of the inhomogeneous derangements 

(6.5) P(n;5,7?)= ^ (J + z)d'-°PW(5 - ^)<=-<-),y^y^W. 

7reI5(Ti) 

Lemma 6.1. For any k,j£ [m] such that k ^ j we have 

(6.6) P+{n;5,7]) - P+{n;5,7]) = 26ink - nj)P{n;5,7]) 

+ rikiuk + + l)Pk (n; S, t]) - nj{nj + i] - 1){5^ + l)Pr (n; 6, rj), 

and the boundary condition ()2.1ip with Xj = 1 for all j and AqCi = Ai. 
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Proof. For j G [m] let n* = (ni, . . . ,nj,l, Uj+i, . . . , Um)- Following the proof of Lemma ISTTl we 
obtain 

(6.7) P{n*;6, rj) = P+{n- 6, 7]) + 2nj6P{n; 6, rj) + (5^ + l)nj{nj + r] - l)Pr (n; 5, r?). 
Subtracting the last equation from (j6.7p with j = k yields (j6.6p . □ 

By the method of separation of variables we can solve (|6.6p and obtain the following result. 
Theorem 6.2. We have 

(6.8) P{n; 6, r/) = / TT (x)u;(x)dx. 

j^w{x)dx Jj^j-J^ 

Proof. Clearly (j2.8p reduces to (j6.6p when Aj = 1 for all j, and 

An = l, Bn = -{5 + 2n)r], = n(r/ + n - 1)(1 + 5^) for all n > 0. 
From Lemma l6. II and Theorem 12.31 we deduce (16. Sp . □ 



This formula was first given by Zeng [35], and later generalized by Kim and Zeng 

7. Linearization coefficients of ^-Hermite polynomials 

The continuous g-Hermite polynomials Hn{x\ q) are generated by 

(7.1) Ho{x\q):=l, Hi{x\q)=2x, 2xHn{x\ q) = Hn+i{x\ q) + {1 - q'')Hn-i{x\ q), n > 0, 
and have the orthogonal relation 

(7.2) / Hn{cos 9\ g)i/m(cos 6*1 q)v{cos e\ q) d9 = (q; q)nSmn, 
Jo 

where 

.(cos0|g) = ^^(e2^^e-2^^g)oo. 



If we rescale the g-Hermite polynomials by 

,1 

'2 



Hn{x\q) = Hn{^ax\q)/a^, a= y^l- q, 



then ()7.ip reads 

xHn{x\ q) = Hn+i{x\ q) + [n]qHn-i{x\ q), 
and the orthogonality relation (|7.2p becomes 

(7.3) / i?n(x|g)i7r„(x|g)'y(x|g)(ix = n!g(5„ 

Here n\q = ((?;(?)n/(l — 9)"" and 



oo 



(7.4) ^(x| q) = J^li + (2 _ ^2(1 _ ^))^. + g2.}_ 

V^l - (1 - g)x2/447r ^^^^ 
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Given a perfect matching M (or more generally, a set partition), a pair of arcs (61,62) of 
M is said to cross if 61 = (i, j), 62 = {k,l), and i < k < j < i. The number of arc crossings 
in M is denoted by cr(M). For instance, if M is the matching drawn in Figured! we have 
cr(M) = 5. Let 

cr{M) 



(7.5) 



K{n\ 



Me/C(n) 



For any nonnegative integer n we set 



(7.6) 



1 



n. 



1 



l + q + ---+q 



n-l 



Lemma 7.1. For k,j E [m] and k ^ j the polynomials K{n\ q) satisfy 
(7.7) K+{n\ q) - K+{n\ q) = [nk]gK~{n\ q) - [n,],Kr{n\ q), 

and the boundary condition ()2.1ip with Xj = 1 for all j and AqCi = Ai. 

Proof. Let u = ui + ■ ■ ■ + uj + 1. The matchings in /C^(n) (resp. /C^^^(n)) can be divided 
into two categories: 

• the integer u E Sj (resp, u E Sj^i) is matched with the ith element u + i in Sj^i 
(resp. u — i in Sj), from left (resp., right), with £ E (resp. i E [uj]), then the 

corresponding arc crosses each of the £—1 arcs of which one vertex is u + 1 (resp. u — t) 
with 1 < t < £ — 1. An illustration is given in Figure [3|^a) (resp.. Figure [3|^b)). Hence 
the generating function of such matchings is 

{l + q + --- + q''^+^~^)Kr^^{n\q) (resp. {1 + q + ■ ■ ■ + q'^^~^)K- {n\ q)); 



1 crossings 




1 crossings 



u u + 1 



u + t 




(b) the blocks Sj and Sj+i in IC^,-^{n) 



(a) the blocks Sj and Sj+i in ICJ {n) 

Figure 3. Crossings in an inhomogeneous perfect matching 



• the integer u is matched with an element not in S'jUS'j+i, let i?ti(n| q) be the generating 
polynomial of such matchings. 

It follows that A'+(n|g) = [njj^i]qK~_^^{n\q) + Ru{n\q) and K'^j^^{n\q) = [nj\qK~{n\q) + 
Ru{n\ q). By subtraction we obtain (j7.7p for adjacent k and j. The general case follows from 
the simple identity Uk — Uj = X]j=j (i*i+i — Ui) for any integers j and k such that j < k. □ 

Theorem 7.2. We have 



(7.8) 



K{n\ 



~ m 

/ v{x\q)T] HnAx\q)dx. 
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Proof. Clearly ()2.8I) reduces to ()7.7p when Aj = 1 for all j, Bk = 0, Cfc = [k]q for all A;, and Ak 
is a constant independent of k. From Lemma l7. II and Theorem 12.31 we deduce (j7.8p . □ 



Remark 7.3. The representation (j7.8p is due to Ismail, Stanton and Viennot [18]. Three 
different proofs were later given in [21 Hm I26j . ^ds we can see, the new proof of (|7.8p given 
above parallels our proof in the case q = 1- 

Note that K{n\Qi) is the number of perfect inhomogeneous matchings of [n] without crossings 
and Hn{x\Qi) is the n-th Chebyshev polynomial of the second kind Un{x). Hence, letting g = 
in Theorem 17.21 we obtain the following result, due to de Sainte-Catherine and Viennot [8]. 

Corollary 7.4. The number of perfect inhomogeneous matchings of [n] without crossings is 
given by 

(7.9) K{n\{)) = - C UnAx)---Un^{x){l-x^f/^dx. 

Another generalization of the above corollary was given by Kim and Zeng |25j . 

8. Linearization coefficients of q-Charlier and (7-Laguerre polynomials 

8.1. Al— Salam-Chihara polynomials. Since our g-Charlier and (^-Laguerre polynomials are 
two rescaled special Al-Salam~Chihara polynomials, we first recall the definition of these 
polynomials. The Al-Salam-Chihara polynomials Qn{x) ■= Qnix]ti,t2\ q) may be defined 
by the recurrence relation [27\ Chapter 3]: 



■1) 



\Qo{x) = l, Q_i(x) = 0, 

[Qn+l{x) = (2X - (ti + t2)q'')Qn{x) " (1 " g")(l " tli2g"-')Qn-l (x) , n > 0. 

Let Qnix) = 2"'pn{x) then 

(8.2) XPn{x) = Pn+l{x) + \itl+ t2)qyn{x) + ^(1 " 9")(1 " tit2q'^^^)Pn-l{x) ■ 

They also have the following explicit expressions: 

^ . . .1 ^ {tit2; q)n , ( q-"^ ,tiu,tiu-^ 

Qn{x; ti,t2\ q) = 3<?'2 



n V il*2,0 



/, . -n , ( q '\t2U ^ 1 



{t2U ■ q)nU 2^1 ( ^-l^-„+l^ 



q;t^^qu ^ 



where x = ^^-^ — or x = cos 6 if u = e*' 



The Al-Salam-Chihara polynomials have the following generating function 

r^f4. \ ^ + + I \ {tlt,t2t; q)oo 

G{t,x) = 2^Qn{x;ti,t2\q)— - 
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They are orthogonal with respect to the hnear functional L: 

(8.3) = — r cos" Q , ^5'^'^''f'''''r''' ^1°" , dQ, 

where x = cos 9. Equivalently, the Al~Salam-Chihara polynomials Qn{x',ti,t2\ q) are orthogo- 
nal on [—1, 1] with respect to the probability measure 

{q,tit2;qU ^ 1- 2(2x2 



(8.4) W'^i''2''^^°° -Q 



^J-Q [1 - 2xtiq>' + tfg2fc][i _ 2xt2q'' + t^q^''] VT^' 

As in [21 [26], we shall consider the g-Charlier polynomials C„(x| q) := C„(x, a, b, c\ q) defined 
recursively by 

(8.5) Cn+i{x\ q) = {x-c- h[n]q) Cn{x\ q) - a[n]qCn-i{x\ q), 

where C_i(x| q) = and Co{x\ q) = 1. Comparing with ()8.ip we see that this is a rescaled 
version of the Al-Salam-Chihara polynomials: 

(8.6) a(x| q) = (-r^Y^'QnCn\f^{x - c - -^); , .01 q] 

\l-qj V2V a 1-q' y/a{l - q) ^ 

We define ui{x) and vi{x) by 

, . _ 1^ 2 _ c{l-q) + b 62 ^ c2(l - g)2 + 2(1 - g)(6c - a) 



2a a 2a(l — q) 

8.7 ^ ^ ^' 

^U^) = o\/— [x-c 



2 V a V 1 
The moment functional for C„(x| q) is 



(8.8) C^{f)^^'''^-' I'-' 



27r 2 V a 

^+ — [1 - 2'ui(x)g^ + g^'']/(2;) dx 



X 



n 1 + 2vi{x)q^/{y/a{l-q)) + (72^/0(1 - (?) y/l-viix^ 



where 



b 

A± = c-\ ± 2 



As in [21], we shall consider the (7-Laguerre polynomials I/n(2;| 9) := Ln{x,y\ q) defined by 
the recurrence: 

(8.9) Ln+i{x\ q) = {x- y[n + l]g - [n]g)Ln{x\ q) - y[n]gLn-i{x\ q), 

with the initial condition L-i{x\q) = and Lo(x| (?) = 1. Hence these are the re-scaled 
Al-Salam-Chihara polynomials: 
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One deduces then the exphcit formula: 

k-l 



fc=0 -I \- \q j=o 



Define U2{x) and ^2(2;) by 
^8 19^ ( \ 2 {l-q){l + y) , j/' + l . . Q-1 y + 1 

^'•''^ ^^(^) = ^^^ ^ ^ + ^^(^)= v^^+ v^' 

Then the moment functional in this case is 
{q,q;q)oo l-q 



i.l3) C2U) 



2vr 2^ 

^+ -pT- [l-2n2(x)g^ + g^^]/(x) 

B_ [1 - 2V2 + ^^Vy] [1 - '^V2 {x)qk+^ ^ + q^^+^y] 

(1 ± ^/y)' 



where 

(8.14) S± 



1 



For the combinatorial approach to the linearization coefficients, the g-Hermite and g-Charlier 
cases were proved by first combining the combinatorial models for the polynomials and mo- 
ments to obtain a messy sum, and then using a killing involution to reduce it to some nicer 
models, [811181 [26]. However, this approach seems difficult to deal with the g-Laguerre case. So, 
a recursive approach based on the symmetry is used in [2Tj, but such a proof for the (7-Charlier 
polynomials is new. 



8.2. Linearization coefficients of g-Charlier polynomials. Recall that if vr is a partition 
of [n], an arc crossing of vr is a pair of arcs (61,62) such that 61 = (i,j), 62 = (fc,^), and 
i < k < j < £. For instance, if tt is the partition drawn in Figured] (resp., in Figure [2|), then 
cr(7r) = 2 (resp., cr(7r) = 6). We let cr(7r) denote the number of arc crossings in vr. 

For each partition vr G n„ we define the weight 

(8.15) W{7r) = abl(7r)^tr(7r)^sg(^)^cr(7r)^ 

where bl(7r), sg(7r) and tr(7r) are respectively the numbers of blocks, singletons and transients 
of TT. Here, a singleton is just a block of size 1 and a transient is an element which is neither 
the least nor the greatest element in a block of vr. 

Consider the enumerative polynomial of inhomogeneous partitions 

(8.16) T{n\q) := J'{n;a,b,c\q) = ^ 'w{tt). 
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Note that by the general theory of orthogonal polynomials, the three-term recurrence re- 
lation (I8.5P and Proposition 4.1 in [22] imply that the linear functional L\ has the following 
combinatorial interpretation: 

(8.17) C^ix-) = J2 

7ren„ 

To find the partial difference equations satisfied by J-{n\ q) we need the following key result. 
Lemma 8.1. The polynomials J-'{n\ q) are symmetric with respect to the permutation of indices 

We postpone the proof of this crucial lemma to Section [TTl 
Lemma 8.2. For j G [m], the polynomials J-{n\ q) satisfy 

(8.18) q) = F{n*\ q) - b[nj\qF{n\ q) - a[nj]gTj {n\ q), 
where n* = (1, rii, . . . , n^). 

Proof. By Lemma |8. 11 we can suppose that j = 1. Hence, it suffices to check that 

(8.19) w{7r)=T+{n\q) + b[ni]gT{n\q)+a[ni]gT{{n\q). 

where w{tt) = a'^l('^) since sg(7r) = for any vr G V{n). 

Given a partition vr E V{n*), we denote by ri the integer i > 1 which is connected to 1 by 
an arc. We classify the partitions in V{n*) into three categories according to the value of ri 
(The reader is suggested to draw diagrams as we do in the proof of Lemma 17. 1|) : 

(a) ri > ni + l\ such partitions are exactly the partitions in 'Pi{n), whence the enumerative 
polynomial of such partitions is J-^{'n\ q). 

(b) 2 < ri < ni -|- 1; then the arc (1,ti) crosses with each of the ri — 2 arcs of which 
one vertex is t with 2 < £ < ri — 1. Suppose {l,ri} is a block of tt (resp., is not 
a block of vr). Summing over all ri = 2, 3, . . . , ni -|- 1, it is readily seen that the 
enumerative polynomial of such partitions is Yl^l=2 aq^^^'^ [n] q) = a[ni\qTi {n\ q) 
(resp., En=2 hq^^-^Fi{n\ q) = 6[ni],^(n| q)). 

Summing up the above three cases we obtain (I8.19p . □ 

The following result is due to Anshelevich [2] and a combinatorial proof was later given by 
Kim, Stanton and Zeng |26j . 

Theorem 8.3. For m > 1 and ni, . . . , Um > 0, we have 

(8.20) T{n\ q) = Ci (C„, (x| q) ■ ■ ■ ?)) • 

Proof. For j, k € [m] and j ^ k we deduce from (j8.18p that 

(8.21) -F/(n|g) - J"+(n|g) = {[nk]g - [nj]g)T{n\q) - a[nj]gTj' {n\q) + a[nk]gJ'- {n\q), 

and the boundary condition (j2.1ip with Xj = 1 for all j and AqCi = Ai. The result then 
follows by applying Theorem 12. 3[ □ 
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Remark 8.4. When q = 0, the polynomials C„(x| 0) are the so-called perturbed Chebyshev 
polynomials of the second kind and J-{n\ 0) is the enumerative polynomial of inhomogeneous 
partitions of [n] without any arc crossings. 

Remark 8.5. In view of Lemmas \8.1i \8.^ and Theorem \2.6\ we can also prove the above 
theorem by checking (j8.20p for the special n = 1™ := (1, . . . , 1). As Ci{x; q) = x — c, the latter 
identity reads 

(8.22) ^(I'^lg) =£i((x-c)"^) (m>l). 
By the binomial formula, this is equivalent to 

(8.23) £i(x-) = ^ r7 c^^^(l™-^|cz). 

A:=0 ^ 

In view of the combinatorial interpretation of the moments ()8.17p and the definition (I8.16P 
the latter identity is obvious if we enumerate the partitions tt of [m] by the weight (18.150 and 
according to the number of singletons. 



8.3. Linearization coefficients of (7-Laguerre polynomials. For a € &n the number of 
crossings of a is defined by 

n n 

(8.24) cr(a) = < i < a{j) < a{i)} + #{j\j > i > a{j) > a{i)}. 

i=l i=l 

Note that the hnear functional £2 has the following combinatorial interpretation [21j : 

(8.25) £2(x") = 

o-66„ 

Consider the enumerative polynomial of inhomogeneous derangements 

(8.26) I{n\q):=I{n;y\q)= yexc(.)^crM_ 

Lemma 8.6. The polynomials I{n;y\ q) satisfy 

(8.27) //(n| q) - /+(n| q) = {yq + 1) (K], - [n,],) I{n\ q) 

- y[nj?qlj{n\ q) + y[nk]llk {n\ q), 
and the boundary condition (j2.1ip with Xj = 1 for all j and AqCi = Ai. 

Proof. It is proved in [211 eq. (38)] that 

(8.28) I+{n\ q) = I{n*\q) - {yq + l)[n,],/(n| q) - y[n,]llr{n\ q), 

where n* = (1, ni, . . . , n^). Replacing j by k in the above equation and then subtracting the 
resulting equation from the above one we get (|8.27p . The boundary condition is obvious. □ 

The following result is due to Kasraoui, Stanton and Zeng |21j . 
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Theorem 8.7. We have 

(8.29) I{n\q)=C2{LnMQ)---Ln„Xx\q))- 

Proof. Clearly (j2.8p reduces to (|8.27p when Xj = 1 for all j, and 

An = l, Bn = -iy[n + l]g + [n]g), C„ = gH^, n > 0. 

From Lemma 18.61 and Theorem 12.31 we deduce (j8.29p . □ 

Remark 8.8. In the above proof, we do not require the combinatorial interpretation of the 
moments (I8.25p . which was needed in [2T] . 

9. More integrals of orthogonal polynomials 

In this section, for a sequence of orthogonal polynomials {pnix)}, we shall consider integrals 
of type 




no e N, 



where fi is an orthogonal measure for {pnix)}- 



One important tool used in this work is MacMahon's Master theorem, \29\ Vol.1, pp. 93-98] 
and its /3-extension due to Foata-Zeilberger |12j . which we now recall. 

Let Vm be the determinant det(5jj — XiOij) (1 < i,j < m). The MacMahon master theorem 
asserts that the coefficient of the expansion of ^ is equal to the coefficient 

of x'|'^X2^ • • • x^™ in 

m 

(9.3) n ( 

k=l 

It will be convenient to restate this in a slightly different form. Let C(m) be the set of 
rearrangements of the word 1"^ . . . m""" . For any rearrangement 

7 = 7(1, 1) • • •7(l,"'i) • ■■li'm, 1) . ..j{m,nm) G C{m), 
we associate the weight 

Then, the coefficient of x"^ ^ • • • in (|9.3p is equal to the sum of all the ^^(7) with 7 
running over all the elements in C{m). On the other hand, each sequence n = (ni, . . . , n^) of 
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positive integers defines a unique mapping x from [n] to [m] given by x{j) = ^ if 3 £ Si. For 
each permutation vr G (3(n) we let 



^tt) = n"x(i),x(^(i))- 

Clearly, to each rearrangement 7 in C(m), there corresponds exactly ni! • • • rim! permutations 
vr in ©(n) with the property that 'w{t:) = w{'y). Therefore, the coefficient of x^^x^ 
(j9.3p is also equal to 

7re6(Ti) 

The MacMahon Master theorem can now be restated as 



m^na ...^„„ in 



X 



ni 



1 . . . 



,1 



>0 " 7re6{Ti) 

The /3-extension of the MacMahon Master theorem [12] reads as follows. 
Theorem 9.1. We have 



(9.4) 



E 



'r-"l . . . 'r'' 



nil - ■ ■ n-m^. 

ni,...,7im>0 7re6(n) 



Now, we consider the determinant 



A. 



m.+l •" 



1 

-X2 

~Xm 
-Xq 



-CXl 

1 



-xo 



-CXl —CXl 

-CX2 -CX2 

1 CX Yfi. 

-Xq 1-Xq 



The proof of the following determinant formula is left to the reader. 
Lemma 9.2. Let a and b be any variables in a commutative ring. Then 



xi a 

b X2 



b b 
b b 



a 
a 



a 
a 



Xn-l a 

b Xr, 



a(pn{b) - b(f)n{a) 



where 4>n{x) = (xi — x)(x2 — x) 
limit (/)„(a) n + a YJj=i 7^ 



[xn — x) . When a = b, the right side should be taken as the 



(9.5) 



Applying the above lemma to A^+i we obtain 
1 



A, 



m+l 



1 



[c(l + Xi) • • • (1 + Xm) - (1 + CXl) •••(! + CXm){l - (1 - c)xq)] 
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Therefore, denoting the elementary symmetric functions of the indeterminates xi, . . . ,Xm by 
ei{x), . . . , em{x), we have 



(9.6) 

Let 
(9.7) 



A^+i = 1 - '^{c + ■■■ + c^-^)ek{x) - xo JJ(1 + 

k=2 3=1 



^H(no,n) = (-l)E7=i"^- 



^"0 



r(a + l)fjL^ 



A main result of this section is the following theorem. 
Theorem 9.3. The integrals {^4^"^ (no, n)} have the generating function 

(9.8) ^^"H^o,n 

no,...,nm>0 



X. 







x: 



m 



no! Um'. 

m 



= 1 - xq JJ(1 + Xj) - e2{x) - 2e3{x) (m - l)e„,(£c) 

Moreover, we have the following combinatorial interpretation: 
(9.9) A^''\no,n)= ^ (a + l)^^^^, 

7r66*(Ti) 

where 6* (n) is i/ie sei of permutations of Sq U ■ ■ ■ L) Sm such that all the elements in box j 
should not stay in the original box after permutation for 1 < j < m and the objects in box 
are not restricted. 

Proof. We use (|4.3p to see that 



no,...,nm>0 



3=0 



— rn(l + ^i)"""' / exp(-2;(l-xo-5^Xfc/(l + Xfc)))x"dx 

1=1 fc=l 



r(a + i) J-J- 

m 



1 



fc=i 



+ Xk 



which reduces to the right-hand side of (j9.8p after some simplification using the following 
identity, which was proved in [6], see also [H (2.8)], 

(9.10) J](l + tj) [l - X] YTf] " ^ " ""'^"^^ - 2e3(a^) (m - l)en,{x). 

j=i j=i + J 

This proves ()9.8p . The rest of Theorem 19.31 follows from the /3-MacMahon Master theorem and 
M. □ 
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Remark 9.4. When a = 0, A^^\nQ,n)/no]nil ■ ■ ■ rim^. can be simply interpreted as follows: 
we have boxes of sizes no,ni,...,nm and box j contains nj indistinguishable elements and 
we arrange the contents such that no object in box j stays in its original box when 1 < 
j < m with no restriction on box number 0. The number of possible rearrangements is 
A(°)(no,n)/no!ni! • • -n^!. 

Corollary 9.5. We have 

(9.11) A^^>{m,n,s) = m\n\s\2_^' 

j>o 

Proof. By (j9.3p we have the generating function 

y A(")(m,n,s)^^^ = y^ 

^-^ m\ n\ S\ \1 — Xn — XiX2 — XiXn — X2Xn — XiX2Xn\°'+^ 

m,n,s>0 

Since 

1 —xi —xi 



n + j — mj \ m 



V 



-X2 1 -X2 
-Xq -Xq I- Xq 



1 — Xq — X1X2 — XiXq — X2X0 — X1X2XQ 



by the MacMahon Master theorem, Theorem 19 -H we see that A^^\m,n, s) is given by the 
coefficient of x^Xj Xq in (x2 + xo)"^(xi + xo)"(xi + X2 + xq)*, which is equal to the claimed 
expression. □ 



Motivated by the numbers A^'^\no,n) we consider the following generalized linearization 
coefficients of Meixner polynomials: 



(9.12) B^^\no,n) = (-1)^^=1 "^■c-"°(l - c)^+"° 



^^(^_l)...(^_no + l)^^n^-.(^;/5'^)- 



x=Q j=l 

Theorem 9.6. The integrals {B^^^uq, . . . ,nm)} have the generating function 



(9.13) B(^Hno,n)ll 



no,...,nm>0 



3=0 



ml 



1 - Yl c(l-c-i) """^"^^ - xollil + Xj/C) 



1-k 



k=2 



:(l-c-i) 



Moreover, we have the following combinatorial interpretation: 
(9.14) B^^\no,n)= ^ ^cycW^-cxcW^ 



where 6*{n) is the same as in Theorem W. 
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Proof. We use ()5.3p to see that 



no,...,nm>0 



nr=i(i + ^j) - (c + (1 - c)xo) nr=i(i + ^./c)^ 



^(1 + (1 _ c)xo/c)-^l^(l - cf + + ^i) 

x>0 ' j=l 



1 - c 

This gives (j9.13p after simphfication. 

Comparing with (|9.6p we see that the (3 = 1 case of (j9.14p comes from the MacMahon Master 
theorem associated with the matrix (ajj) with an = 0, Oij = 1/c for j > i and aij = 1 for j < i. 
The general case fohows from using the /3-extension of MacMahon's Master theorem. □ 

Remark 9.7. For the Charlier polynomials we have a similar result for the integral 
(9.15) C(")(no, n) := x(x - 1) • • • (x - no + 1)^^ f[ C^^^ix). 

x>0 ' j=l 

A straight computation shows that 

no ni rim 



(9.16) C(''\no,n)^^ 



no,...,n,n>0 

= exp {a[xo + xoei{x) + {xq + l)e2(a;) H 1- (xq + l)em(a;)]) 

VFe apply the exponential formula to see that 

bl{7r) 

7r6'P*{no,Ti) 



(9.17) C7W(no,n)= ^ a 



where V*{no, n) is the set of partitions o/ 5o U 5i U • • • U Sm such that each block is either a 
singleton of an element in So or inhomogeneous, i.e., no two elements of Sj (0 < j < in) can 
be in the block. 

It is clear that Theorem 19.31 is the limit c — > 1^ of Theorem 19.61 Similarly we have the 
following analogue of Corollary 19.51 

Corollary 9.8. We have 

B<..K„.,.»i„i.E(';)(„^;_J (-:-■).--■. 

Corollary 9.9. We have 

(9.18) = [I) + ^)n-fc(-l)'Mfc(x; (3, c). 
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Proof. Let = ^^^q c(n, A;)Mfc(x; /3, c). Using the orthogonality (I5.2|) we obtain 

(1 - ^""M.ix; /3, c)M^c- = c(n, k)^-^. 

Comparing with (|9.12p we see that the left side is equal to (— l)'^c"'(l — c)~"'B^^\k,n). It re- 
mains to compute B^^\k, n), which, by Theorem l9.6l is the coefficient of -^r^ in the expansion 

n>0 n,A:>0 ^ ^ 

Hence B^^\k,n) = jj^c"''. This yields the desired result. □ 

Let if be the linear functional defined by <p{f{x)) = j-^f{x)dii{x). Then the integral (j9.ip 
contains the following four special cases: 

(1) the evaluation of fix"") corresponds to the moments, 

(2) the evaluation of ^^(11^=1 Prij (^)) corresponds to the orthogonality, 

(3) the evaluation of if{x^pk{x)) combined with the orthogonality corresponds to the co- 
efficient Cn,k in the expansion x" = X]fe=o (^n,kPk{x), 

(4) the evaluation of ^{YVjLiP^ji^)) corresponds to the linearization coefficients. 



Since Aqx = pi{x) — Bq, we have 



(^oxr = f;(^)(-Bo)^"' 



Piix), no G N. 



Therefore, 



(9.19) ^((^ox)"°n^'".(^)) =E 7r"^°)""'^(^^^'^)'n?'".(^))- 

i=l 1=0 ^ ^ j=l 

We can deduce the combinatorial interpretations of the integrals ()9.19p for the orthogonal 
Sheffer polynomials and the three g'-analogues from the combinatorial interpretation of the 
corresponding linearization coefficients. 

For example, as Hi[x) = 2x, it follows from Theorem 13.21 that 



(9.20) 2-("o+"i+-+"™)/2 / l^(2x)"« JJi?„^(x)d: 



e '^^ 



vr 



is the number of perfect inhomogeneous matchings in /C(n) with 

n = (1, . . . , l,ni,n2, . . . ,nm)- 



no 
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For the Laguerre polynomials we have x = —L^^\x) + a + 1, so 

We can easily recover the combinatorial interpretation ()9.9p in Theorem 19.31 from the above 
equation and ()4.13p . 

For the Meixner polynomials we have ^—jrX = P — Mi{x; /3, c), so 

(9.22) ^('^)(no,n) = c-"«(l-c)^+"o 

a;=0 ^' j=l 

Using Theorem 15.2^ we see the following combinatorial interpretation 

(9.23) B^^\no,n)= ^ ^cycW^-excW-excoW^ 

where &*{n) is the same as in Theorem 19.31 and exco(vr) is the number of excedances of two 
elements in Sq, i.e., exco(vr) = \{i £ Sq : 7r(z) G 5o and 7r(z) > 



10. Laguerre and Meixner polynomials revisited 

Recall [171 p. 100] that the Hermite polynomials can be viewed as special Laguerre polyno- 
mials since 

H2n+1/2±1/2{X) = (-l)"22"n!(2.)V2±l/24±l/2)(^2), 

Therefore the integral in (j3.3p is a special case of the integral 



(10.1) Wj^k{m; a, /3; m, n) := 



(_l)ELi™i+ELi"'- 
r(a + l) 



r J 



r k 



i=l 



■r=l 



dx, 



where m = (mi, m2, . • • , rrij) and n = (ni, n2, . . . , n^). 

Li this section we study the combinatorics of the integrals of the type in (|10.ip and their 
discrete analogues which result by replacing the Laguerre polynomials by Meixner polynomials. 



SEPARATION OF VARIABLES IN COMBINATORICS 



33 



Theorem 10.1. Let ei{x), i = 0,1, . . . ,j + k, be the ith elementary symmetric polynomial of 
xi, . . . , Xj+k- The integrals {Wj^ki^; a, P; m, n)} have the generating function 



(10.2) Yl Wj,k{m;a,(3;m,n) 

m, nii , nr>0 

k 



Xq X^ X, X- 



X 



■j+k 

ml mil rrijl nil n^l 



r=l 



j+k j+k 

l-xoll{l + x,)-Y,ii-l)ei{ 

1=1 1=2 



H -a-1 



Proof. Apply the generating function (|4.3p to see that the left-hand side of (|10.2p is given by 

dx 



J 1^ r-oo a / _ _ 

(10.3) \{{i + x.)—'X{{i + x,^r)-'-' f7;7Tn^"P -^ + ""° + Eti^ 

i=l r=l Jo L) \ 1=1^^^ 

k 



11(1 + x,r-^ + x,+.)-^-i [i-xo-Y.- 



j+k 



+ Xl 



Xl 



i=l 



r=l 



+ Xl 



-a-1 



This establishes (110. 2p after some simplification using (I9.10p . 



□ 



Corollary 10.2. The numbers {Wj^k{m;oi,(5;m,n)} are positive when a > — 1 and a — j3 is 

a nonnegative integer. 



Assuming that a — /3 is a positive integer N , we can give a combinatorial interpretation for 
Wj^k{fn', Oi, 13; m, n). Let (3^(n) be the set of {k + l)-tuples (vr, /i, . . . , fk) such that 

• cr is an inhomogeneous permutation of U U • • • U S'j U S*_^_i U • • • U S*_^_f,, where 
5*0 ^ 5*0 and S*_^_j. C Sj+r for r = 1, . . . , /c. 

• fr ■ Sj+r \ Sj^^ — )■ [N] is an injection for r = 1, . . . ,k. 



From Theorems 19.31 and 110.11 we deduce the following combinatorial interpretation: 

(10.4) Wj^k{m;a,P;m,n) = ^ (a + l)^^^^. 

{7r,/i,...,/fe)66;,{n) 

Motivated by the numbers Wj^k{m]0(, I3;m,n) we consider the following generalized lineariza- 
tion coefficients of Meixner polynomials: 

(10.5) Yj^k{m;a,P;c;m,n) := (-1)^=1 "-c-™(l - 0)"+"^ 

X ^x(x-l)---(x-m+l)^^ J[M^^{x;a,c) JJ M„,(a:;; /3, c) 



x=0 



i=l 



r=l 



where m = {m\,m2, . ■ ■ , rrij) and n = (ni, 712, . . . , n^). 
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Theorem 10.3. The integrals Yj^k{rn]a, j3;c]m,n) have the generating function 



(10.6) ^^>K"'/5;c;m,n)-^n;^ n 



m! mil ^ rir'- 



I 



m,mi,nr>0 i=l r=l 

k [ ^^'^ I _ . 

n(i + xj+rT-^ 1 - ^ ei{x) - xo n(i + ^^/t 



1=2 ^ ' i=\ 



Proof. Applying (|5.3|) to see that the left-hand side of (|10.6p is 

Yii + ii-c)xo/cr^^ii-cr 

x>0 

j k 

+ xjcni + 11(1+ x,+r/cni + ^ 



X 

i=l r=l 
j k 



j+k 



1 + Xi/c 



. , - + Xi 
1=1 



= (1 - cT n(i + ^0"" 11(1 + 1 - (c + (1 - c)xo) n -i 

i=l r=l '- i=\ 

This establishes (jl0.6p after some simplification using ()9.10p . □ 



In the same vein, assuming that a — /3 is a positive integer N , Theorems 19.61 and 110.31 implv 
the following combinatorial interpretation: 

(10.7) Yj^k{m;a,P;c;m,n) = ^ 

(vr,/i,...,/fe)66-V(ri) 



Q,cyc(7r)^-cxc(7r)_ 



Note that Theorem 110.31 shows that the numbers Yj k{rn; a, /3; c;m,n) are positive when 
a — /? is a nonnegative integer. 

11. Proof of Lemma [871} Symmetry of T{n\q) 

Recall that n = (ni, . . . ,nm) is a sequence of positive integers and n = ?7,i + • • • + Um- 
Clearly we need only to prove the invariance of J-{ti\ q) for the two following permutations of 
the indices n^-'s: the transposition exchanging 1 and 2, and the cyclic permutation mapping i 
to i + 1 (mod m) for i = 1, . . . ,m. Moreover, since sg(7r) = and tr(7r) = n — 2bl(7r) for any 
partition vr G 'P{n), we see that Lemma |8. H is equivalent to the following result. 



Lemma 11.1. We have 

-|bl(7r)^cr(7r) 

TreV^n) n£'P{n2,n3,...,nm,ni) 

-jbl(7r)^cr(7r) 

7rGP(Ti) ne'P{n2,ni,n3...,nm) 



(11.1) Y a'^'^^U''^''^ = E 

7rg7'(n) 7rg7'(n2,n3,... 

(11.2) Y a'^'^^U''^''^ = E 
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For a positive integer k such that k < n, we introduce two sets of inhomogeneous partitions: 
«P„ ■.= V{k,l^_^), P« ■.= Vil^_^,k). 

n—k n—k 

In other words, a partition vr of [n] is in ^^^Vn (resp., Vn^) if and only if it has no singleton 
and there is no arc in vr joining two elements in [1, k] (resp., [n — k + ^,n]). For instance, the 

two partitions vri and tt2 drawn at the top of Figure [Hare in ^^^T^ia and 'P['^K We first show 
that the following result implies (jll.ip . 

Proposition 11.2. For any positive integer k, there is a bijection : ^'^^Vn ^ such 
that for any vr G '•^^"Pn? we have 

(I) for k < i < j , the pair {i,j) is an arc of vr if and only if the pair (i — k,j — k) is an 

arc o/^>n,fc(vr); 
(II) bl($„,fc(7r)) = bl(vr) and cr(^>„_fe(7r)) = cr(vr). 

Indeed, assuming the existence of such a bijection ^n,k with k = ni, as r{n) C the 
property (I) implies that $„,^„j('P(n)) C ^(722,^3, . . . ,nm,ni). Since the cardinality of V{n) 
is invariant by permutations of the rij's and ^n,ni is bijective, we deduce that 

$n,ni CPin)) = P(n2, ns, . . . , Um, ni), 

and then (jll.ip by applying the property (II). 

We now turn our attention to (jll.2p . Define the set of inhomogeneous partitions 

-p(m,n,) :=r{m,n2,l^_^). 

n—n\—n2 

In other words, a partition vr of [n] is in -p^"!'"^) if and only if it has no singleton and there 
is no arc connecting two integers in [l,ni] or in \n\ + l,ni + 712]. For instance, the partitions 

TTi and 7r2 drawn in Figure [6] are, respectively, in vf^^^ and V^f"^ . Similarly, we deduce (|11.2p 
from the following result. 

Proposition 11.3. There is a bijection G^T^'"^^ : -p^'^i.'^z) _^ p^n2,ni) gy.^y^ ^^^^ ^^^y 
TT £ Hn .we have 



(I) for N2 < i < j, the pair is an arc of vr if and only if the pair {i,j) is an arc 

of @^n^'"'^\Tr), where N2 := ni + 712; 
(II) bl(Gi"^'"'^(^)) = bl(^) and cr(ei"^'"'^(^)) = cr(7r). 



Indeed, since P(n) C vt'''"^\ the property (I) of G^T''"'^ implies that Gi"''"'^(P(n)) C 
V{n2,ni,n3, . . . ,nm)- This, combined with the fact that the cardinality of 'P{n) is invariant 
by permutations of the n^'s and Gn""^'""^^ is bijective, implies that 

e(ni,n2)(^(^)) = r{n2,nun3, n^). 

Equation (jll.2p then follows by applying the property (II) of Q^^''^'^\ 

The next two subsections are dedicated to the proof of Propositions 111.21 and lll.3[ 
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11.1. Construction of the bijection ^n,k- Given a partition vr G n„, an element i G [n] 
is said to be minimal (resp., maximal) if i is the least (resp., largest) element of a block 
of vr. The set of the minimal (resp., maximal) elements in vr will be denoted min(7r) (resp., 
max(7r)). For example, for vr = 146/2/3 5, min(7r) = {1,2,3} and max(vr) = {2,5,6}. Note 
that min(vr) n max(vr) = sing(7r) where sing(7r) is for the set of singletons of vr. Let S be a 
subset of X. The restriction of a partition vr = {Bi, B2, ■ ■ ■ , B^} of X on S is the partition 

{Bi n 5, ^2 n 5, . . . , n S} of s. 

The key idea for the definition of the mapping $ri,fc is some appropriate decomposition of 
partitions in 

^'''>Vn and vi^K Let (resp., A^f^) be the set of 3-tuples {T,R,a) where 

• r G Un-k and a G &k, 

• sing(r) C i? C min(r) (resp., sing(r) C C max(r)) and \R\ = k. 

For instance, in Figure [U we have (ri,0,cri) G ^^^^13 and (r2,C, 172) G ^13^- 

We first define two simpler mappings Fn^k '■ 'Prt^ ^rf^ and Gn,k '■ ^'^^'Pn — ^ ^^'^^n- 

• For vr G vit\ set -F„,,fc(vr) = (r, C, fj), where 

— r is the restriction of vr on [n — A;] ; 

— C is the set of elements in vr which are connected to an element > n — by an arc; 

— By definition of vit\ we have |C| = k. Suppose C = {ci < C2 < • • • < c^}, then a 
is the unique permutation in such that {ci,n — k + cr(l)), (c2, n — k + cr{2)),. . . , 
{ck,n — k + cr{k)) are arcs of vr. 

• For vr G ^'^^Vn, set G„^/fc(vr) = {T,0,a), where 

— r G Hn-k is the partition obtained by subtracting k from each element in the 
restriction of vr on [fc + 1, n]; 

— Let M be the set of elements in vr which are connected to an element j < k hy 
an arc. By definition of vit\ we have \M\ = k. Suppose M = {mi < m2 < 
■ ■ ■ < nik}, then O is obtained by subtracting k from each element of M, i.e., 
O = {mi — k, 7112 — k, . . . , ruk — k}; 

— cr is the unique permutation in &k such that (cr(l),mi), ((T(2),m2),. . ., {a{k),mk) 
are arcs of vr. 

The mappings Fn^k and Gn,k are illustrated in Figure [H 

Definition 11.4. Let n be a partition of a set S consisting of positive integers. The depth of 
an element i in it, denoted dpj(vr), is the number of arcs {a,b) in it satisfying a < i < b. 

Definition 11.5. Let a = (t(1) cr(2) • • • o"(n) be a permutation of [n]. A pair 1 < « < j < 

n, is said to be a non-inversion in a if cr{i) < a{j). The number of non-inversions in a will 
be denoted nmy{a) . 

Some useful properties of Fn^k and Gn,k are summarized in the following result. 
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= {1,2,3,5} — ^C = {2,4,5,9} 

ai = 3142 ^^'"^ 0-2 = 3142 

Figure 4. The mappings ^n,k, Fn,k, Gn,k and ^n,fe 

Proposition 11.6. The mappings Fn^k '■ "Pn^ ~^ ^i^^ o^i^-d Gn^k '■ ^^^'Pn -^n O'f^ bijections. 
Moreover, for any tt G vi^\ if Fn,k{'^) = {T,C,a), then 

(11.3) bl(7r) = bl(r) and cr(7r) = cr(T) + ninv(cT) +^^dpj(r), 

iec 

and, for any it G ^^^Vn, if Gn^ki^^) = {T,0,a), then 

(11.4) bl(7r) = bl(r) and cr(7r) = ctc{t) + muv{a) + ^dpj(r). 

Proof. It is easy to see that Fn^k (resp., Gn^k) is a bijection by constructing its inverse (use 
Figure . 

Let S" be a finite subset of positive integers. Clearly, if vr is a partition of S, then each 
block i? of TT is represented by \B\ — 1 arcs. This easily leads to the following result. 

Fact 11.7. The number of blocks of a partition tt of S is equal to \S\ — {number of arcs in tt). 

The first equation in (jll.3p and ()11.4p is just a consequence of the above fact. We now 

turn our attention to the second equation in (I11.3P and (lll.4p . Let tt G vit^ . Clearly, the arc 
crossings in the partition tt can be divided into three classes -Ri(vr), R2{tt) and Rs^tt) illustrated 
in Table [TJ 

They are defined formally as follows: 

-RiW = {(n,ii)(«2, J2) G vr I 1 < ii < ^2 < ii < i2 < - A;}, 
-R2(vr) = {(n,ii)(i2, J2) €TT\l<ii<i2<n-k<ji< J2}, 
Rsi-^) = {(n,ii)(«2, J2) €TT\l<ii<i2<ji<n-k< J2}; 

and satisfy cr(7r) = |i?i(7r)| + |i?2(7r)| + |i?3(7r)|. Suppose -Fn,fc(vr) = (r, C, o"). Then it is easily 
checked that |-Ri(7r)| = cr(T), |i?2(7r)| = ninv((7) and |i?3(7r)| = J2i£C '^Pi^'^) (^^^ Figure H]). 
This proves the second equation in (jll.3p . 
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2, 






1 






1 k n 


1 n - k n 


2 






1 k n 


1 n — k n 


3 






Ik n 


1 n — k 11 



Table 1. Sketch of crossings in Lj(7r) and Ri{'K). 



Similarly, let vr € Vn ■ The arc crossings of the partition vr can be divided into three parts 
Li(7r), L2{n) and L3{n) illustrated in Tabled! defined formally as follows: 

Li{Tr) = {(n, ji)(«2,i2) G vr I < ii < i2 < ii < i2 < n}, 
L2{'^) = {(n,Ji)(«2,i2) G vr I 1 < Zi < Z2 < ^ < ii < J2 < ?^}, 
-^3(7r) = {(n, ji)(«2,i2) G vr I 1 < zi < A: < Z2 < ii < J2 < ?^}, 

and such that cr(7r) = |Li(7r)| + |L2(7r)| + |L3(7r)|. Suppose G„_fc(7r) = (r, 0,cj). Then it is easily 
checked that |Li(7r)| = cr(r), |L2(vr)| = ninv(cj) and |L3(7r)| = XlieO ^Pj(''") (^^^ Figure H]). 
This proves the second equation in (jll.4p . □ 

In view of Proposition 111.6^ to prove Proposition 111.21 it suffices to prove the following 
result. 

Proposition 11.8. For any partition vr, there is a bijection : min(7r) i— >■ max(7r) such that 
dPj(vr) = dp^(j)(7r) for each i G min(7r) and = j for j G sing(7r). 

Proof. It is worth noting that such a bijection was already described in the literature (e.g., 
see Remark 7.2 in |23]). For reader's convenience we recall the construction of V'tt- The 
mapping ipT^ can be nicely illustrated using Motzkin paths. Recall that a Motzkin path of 
length n is a lattice path in the plane of integer lattice from (0,0) to (n,0), consisting of 
NE-steps (1, 1), E-steps (1,0) and SE-steps (1, —1), which never passes below the x-axis. The 
usual way to associate a set partition to a Motzkin path works as follows: to a partition tt of 
[n] we associate the Motzkin path M of length n whose i-th step is NE if z G min(7r) \ sing(7r), 
SE if i G max(7r) \ sing(7r) and E otherwise. An illustration of this correspondence is given in 
Figure [5l 

A basic property of the above correspondence is the following fact |22j . 

Fact 11.9. Suppose M is the Motzkin path associated to a partition vr and let hi he the height 
of the i-th step of M , i.e., the ordinate of its originate point. Then, dpj(7r) = hi if the i-th 
step of M is NE and dpj(7r) = hi — 1 if the i-th step of M is SE. 
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Figure 5. The Motzkin path associated to the partition vr = 
1415/2 3/5 6/71013/8/911/1214 and the mapping Vtt 

We can now describe the mapping ip^. We first set V'ttO) = j for j G sing(7r). Suppose 
©(tt) := min(7r) \ sing(7r) = {oi < 02 < • • • < Or}, C{ir) := max(7r) \ sing(7r) = {ci < C2 < • • • < 
Cr} and let M be the Motzkin path associated to vr. Note that the NE (resp., SE) steps in M 
are exactly the steps indexed by 0{tt) (resp., C{-k)). We then pair the NE-steps with SE-steps 
in M two by two in the following way. Suppose the i-th NE-step (i.e., the Oj-th step) of M 
is at height h. Then, if the first SE-step to its right at height h + 1 is the j-th SE step (i.e., 
the Cj-th step) in M, then we set il^n{oi) = cj. An illustration is given in Figure [5j From the 
construction of "077 and Fact 111.91 it is easy to see that V'tt is the desired bijection. □ 

For {T,0,a) E ^'^^An, we set ^n,k{''',0,a) := {T,tl^r{0),a). Clearly ^n.k is a mapping from 

^'''^An to An illustration is given in Figure HI From Proposition 111.81 we immediately 

deduce the following result. 

Proposition 11.10. The mapping : ^'^^^n — ^ A^^ is a bijection. Moreover, if {T,0,a) G 
^^^An and ^n,k{'T,0,a) = {T,C,a), then we have 

ieC i&O 

Finally, we define the mapping $„,^fc : ^^^P„ — vit'^ by 
(11-5) $n,fc := F'l o o Gn,k- 

This mapping is illustrated in Figure [H Combining Propositions 111.6) and [Tl.lOl we conclude 
that the mapping ^n,k satisfies the requirements of Proposition 111.21 

11.2. Construction of the bijection Q^^'^'^'^ . The key idea for the definition of the mapping 
0(rii,n2) -g gQj^g appropriate decomposition of partitions in 'piC^'"'^^ We first introduce some 
further definitions. For any set K, let Il{K) be the set of partitions of K. 

Definition 11.11. For two positive integers r, s, we denote by V*{r, s) the set of all partitions 
vr of [r + s] such that there is no arc in vr connecting two integers in [l,r] or in [r + l,r + s] 
but vr can have singletons. Thus, we have V{r,s) C V*{r,s). 

Definition 11.12. Let A^n^''^^^ be the set of 3-tuples {{t, A), {j, B),a) where 
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• T is a partition in V* {ni,n2) and A is a set satisfying sing(T) C ^4 C max(T); 

• 'J is a partition in n([A'^2 + 1,^^]) o,nd B is a set satisfying sing(7) C C min(7); 

• the sets A and B have the same cardinality. If k = \A\ = \B\, then a is in Gk- 

For instance, in Figure^ we have ((ti, Ai), (71, cri) G A^^^^ and ((t2, ^2), (72, -B2), 0-2) G 

,(4,3) 

For vr G v}r''"'\ we set H^T^'^'X^) = ((r,^), {-i,B),a) where 

• r is the restriction of vr on [1,A'^2] and A is the set of elements < N2 in vr which are 
connected to an element > A'^2 by an arc; 

• 7 is the restriction of vr on + li ^] and B is the set of elements > N2 in vr which are 
connected to an element < A'^2 by an arc; 

• Suppose A = {oi < a2 < • • • < Ofc} and B = {hi < 62 < • • • < ^k}- Then, a is the 
(unique) permutation in (3fc such that (ai, h^^i-^), (02, &ct{2))i • • (ofc, bo-(fc)) are arcs of vr. 




Bi = {9,10,11,13} 
CTi = 2413 




B2 = {9,10,11,13} 
CT2 = 2413 



Figure 6. The mappings ri"^'"^) and 9^"^'"^^ 

Proposition 11.13. The mapping Hn ^' ; 'pj^ ^' — > An ^' is a bijection. Moreover, for 

anyTreVj{''''^'\ 2/ifi''''"'^(7r) = ((r, A), (7, 5), cj) andk=\A\, then 

(i) bl(^) = bl(r) + bl(7) - A;, 
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(ii) cr(7r) = cr(r) + 01(7) + nmv(o-) + Y^i^^j, dpi(r) + Y.i<^B dPi(7)- 

Proof. It is easy to see that H^^'^^^ establishes a bijection from 'p^"!'"^) ^^"^'"^^ by con- 
structing its inverse (use Figure [6]), and Property (i) is a direct consequence of Fact lll.7[ 
Let vr G ■p^"^'"^^ The arc crossings of the partition tt can be divided into five parts Ci('/r), 
1 < i < 5, iUustrated in Table [21 They are defined formally as follows: 



Ci 




= {(n,ii)(«2,i2) G vr 


1 < n < «2 < ji < i2 < A^2}, 




C2 




= {(n,ii)(«2, ^2) G vr 


N2 <ii <i2 < ji < J2 < n}, 




C3 




= {(n,ii)(«2,i2) G vr 


1 < n < «2 < A^2 < Jl < J2 < 


n} 


C4 




= {(n,ii)(«2,i2) G vr 


1 < < iV2 < i2 < Jl < J2 < 


n} 


C5 




= {(n,ii)(«2,i2) G vr 


1 < n < «2 < Jl < iV2 < J2 < 


n} 



and satisfy cr(7r) 

= ELi l<^i(^)l- Suppose Hr^'^'^'iir) = {{r,A), {7,B),a). It is easily checked 
(use FigureED that |Ci(7r)| = cr(T), |C2(7r)| = cr(7), |C3(7r)| = ninv(a), |C74(^)| = Ei&B^^Pih) 
and |C5(7r)| = X]jg^dpj(r). Altogether, this leads to Property (ii). □ 



C,(7t) 



1 "1 



^2 



N2 



N2 



N2 



1 "1 



N2 



Table 2. Sketchs of crossings in Ci(7r). 



Let 

(11.6) i?(ni, 712) := {(vr, A) : vr G 7^*(?ii, 712) and sing(vr) C A C max(vr)}. 

For instance, the elements (vr, A) and (vr. A') drawn in Figure [7| are, respectively, in -R(4, 6) and 

In view of Proposition 111.13^ to prove Proposition 111.31 it suffices to demonstrate the fol- 
lowing result. 
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Proposition 11.14. There is a bijection ip{ni,n2) '■ -^('^1)^2) — ^ Rin-2,ni) such that, for 
{tt,A) G R{ni,n2), i/ ^(„^,„2)(7r, A) = {7:',A'), then 

(11.7) cr(7rO = cr(^), \A'\ = [A], ^ dp,(7r') = dp,(vr). 

Proof. To any {ir,A) G R{ni,n2) we associate an element {tt',A') in R{n2,ni) as follows: 

• By definition of V*{ni,n2), the arcs of vr are (ii,ip(i)), (^2, jp(2)), • • {ik,3p(k)) for some 
integers A; > 0, 1 < ii < ^2 < • • • < ^fc < f^ii '^i + 1 < ii < ^2 < • • • < ifc < ^2 and some 
permutation p G ©fc. We use z for the complement of i in [1, A'2], i.e., i = N2 + I — i. 
Then, we define vr' as the partition of [l,iV2] which consists of the arcs {jr,ip{r)) for 
1 < r < k. It is clear that vr' G V* {n2,ni). Moreover, we have cr(vr') = nmv{p) and 
cr(vr) = ninv(/9) whence cr(vr') = cr(vr). 

• Since sing(vr) C A C max(vr), we have A = sing(vr) U B with B = {^^(i) < J^(2) < • • • < 
je(^t)} for some increasing sequence (^(s))i<s<t. Suppose / := {«i,«2i • • • ,ik} = {^^i < 
U2 < ■ ■ ■ < Uk}- We then set A' := sing(vr') U B' with B' = {^^(i) < ti£(2) < • ■ • < 'U^(j)}- 
Clearly, we have sing(vr') C A' C max(vr') and \A'\ = \A\. It is also easily checked 
that (in,(,)(vr') = iij,(,j(vr) for s = 1,2,. whence EieB' dPi(vr') = EieBdPi(^)- 
Moreover, since sing(vr') = sing(vr) and dj(vr') = di{7r) for i G sing(7r), we see that 
Eiesing(7r')dPi(7r') = Eiesmg(7r)dPi(vr). Altogether, this implies that T.ieA' ^Pii'^') = 
EieAdPi(7r). 

Set „2)(^i^) = (^'i^O- Then ip[ni,n2) is a well-defined map from R{ni,n2) to R{n2,ni) 
and satisfies ()11.7p . An illustration is given in Figure [71 Besides, it is easy to see that the 
composition V'(n2,ni) ° V'(ni,n2) the identity mapping. This proves that V'(ni,n2) a bijection. 

□ 




3, 5, 6, 7, 8, 10} {1, 4, 6, 7, 8, 9 } = A' 



Figure 7. The mapping '0{ni,n2) 

For ((vr,^),(7,S),<T) G we set 

r^'"^)((vr. A), (7, 5), a) := (^(„„„,)(t, ^), (7, S), a). 

Clearly r^"!'"^) jg mapping from A^^'^^^ to ^4^"^'"^^ An illustration is given in Figure [H 
From Proposition 111.1^ we deduce the following result. 

Proposition 11.15. The 

{{T,A),{^,B),a) G 4"''"^) an^fr^'"^)((vr,A),(7,i3),a) = ((r', A'), (7, S), a), then we have 
cr(r')=cr(r), |A'| = J] dp,(r') = ^ dp,(r). 
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Finally, we define the mapping B^T'''''^ : pi"^'''^) ^ p^"2,ni) 

(11.8) e^r''"^) := O o /^K,rX2), 



An illustration is given in Figure [6l Combining Propositions 111.13] and [11.151 we conclude that 
the mapping 0^"^'"^^ satisfies the requirements of Proposition 111.31 
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